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I. INTRODUCTION 



Interpretation of the data obtained by means of different methods of the positron annihi- 
lation spectroscopy (PAS) is largely based on theoretical understanding of primary physico- 
chemical reactions and processes in the terminal part of the positron track, and on comprehen- 
sion of the nature of the states from which the positrons may annihilate. ^'^ 

Because of extreme rapidity of the primary radiation chemical processes, fundamental knowl- 
edge about them is obtained mainly from theoretical analysis and interpretation of relevant ex- 
perimental data and first of all from pulse radiolysis, permitting to enter the area of pico- and 
(probably, in the near future) femtosecond time scale.^ PAS and especially the experiments in 
external electric fields allow one to obtain important supplementary information about the intra- 
track processes, namely ion-electron recombination, e'^ and e~ thermalization and localization, 
positronium (Ps) formation, influence of various chemical additives and to make conclusions 
about the spatial structure of the e+ track.^ These advances throughout the last two decades 
became possible due to essential instrumental achievements (application of low-energy positron 
beams, better time resolution, accumulation of higher statistics) and theoretical development 
(calculations of cross sections and rate constants of a variety of fundamental processes with 
participation of positrons, investigations of e''"-e~-correlations, slowing down, track effects). 
This progress significantly increased the reliability of PAS. 

First experiments on Ps formation (including those in external electric field) were carried 
out in low-density simple atomic gases, where Ps is formed predominantly according to the Ore 
mechanism. ^'^ It implies that the "hot" positron, e"*"*, pulls out an electron from a molecule M, 
forming Ps atom and leaving behind the radical-cation M"*"': 



This process is the most effective, when the energy W of the positron lays within the interval 
named the Ore gap: 



Here Iq is the first ionization potential of the molecule, W^x is its electronic excitation threshold 
and Ry/2 = 6.8 eV is the Ps binding energy in vacuum. It is believed, that the positron with 
the energy lower than Iq — Ry/'2 cannot pick up an electron from a molecule. When W > Wex 
electronic excitations and ionizations compete with Ps formation and the last one becomes less 
effective. 

Theoretical consideration of the electric field effect on Ps formation in gases was based on 
a solution of the Boltzmann kinetic equation.^ This approach successfully explained the rather 
sharp growth of the Ps yield vs the field in rare gases, molecular hydrogen (H2, D2) and N2 
at moderate densities.^ Electric field hinders energy loss by the positrons and at the same 
time "heats" them up to the energies above the lower boundary of the Ore gap. The simple 
estimation of the field E needed to keep e'^ within the Ore gap can be obtained if we equate a 
gain of the e"*" kinetic energy between two subsequent collisions 



and the average energy loss, (this expression corresponds to the energy loss in elastic 

collisions). Here m is the e"*" mass, M is the mass of a gas molecule and Itr is the positron 



e+* + M^Ps + M+-. 



(1) 



Ig - Ry/2 <W <W, 





(3) 
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transport mean free path (Appendix A). Averaging (. . .)y is carried out over orientations of 
the velocity v of the positron. From such a balance condition one obtains that the average 
e"*" kinetic energy {W)y is approximately equal to [M/mY^'^eEltr- At typical gas densities and 
external fields 10-500 V/cm it falls in the eV region. Thus, rather moderate electric fields are 
able to accelerate positrons up to the energies needed for realization of reaction (1). It is worth 
mentioning that in liquid helium the field dependence of the Ps formation probability, Pps, is 
qualitatively the same as in the gas phase.'' It means that the Ore process contributes to the 
Ps formation in liquid helium. Respective growth of the Pps occurs in fields, approximately 
500 times larger than in the gas phase, which is in accordance with the above estimation. This 
increase of the threshold field is determined by the decrease of Itr, which is inversely proportional 
to the density of helium. 

More detailed theoretical considerations of Ps formation problem in light noble gases require 
a knowledge of energy dependencies of momentum transfer, a^, and annihilation, a a-, cross sec- 
tions of e"*" vs its kinetic energy (see Ref and references therein) . In terms of energy dependent 
cross sections the theory succeeded to explain the most striking experimental observation in 
gases: so-called "shoulder" in the e"*" lifetime annihilation spectra.^ It appears due to existence 
of the Ramsauer minimum in (JmiW) which usually takes place at ~ 1 eV. 

In the 1990's, the Ps formation in noble gases was studied over a wide range of densities in 
the presence of static electric field. One of the motivations for these studies was a significant 
contradiction between available experimental values of Pps in heavier noble gases (Kr, Xe) and 
corresponding predictions of the standard Ore model. Another unanswered question was why 
the application of a strong electric field does not increase Pps up to unity? Probably the origin 
of the effects (as well as a surprising decrease of Pps at higher fields in Xe) is related to the 
formation of locahzed states (due to the self-trapping) by both e"*" and Ps prior to annihilation.^^ 

In molecular gases (H2, N2, CO2, CH4) with low-lying vibrational or rotational levels, Ps 
yield reveals complex behavior with changes in the density, temperature and electric field.^ 
It cannot be convincingly treated in terms of e+-M and Ps-M cross sections. However exper- 
imental data can be fitted by means of combination of the Ore process and recombination 
mechanism. The last one suggests Ps formation in condensed media through combination of 
the thermalized positron with one of track electrons produced by e+ ionization of the medium 
in the terminal part of the e"*" track. ^^^^^ Later on, the Onsager-like formulation of the re- 
combination mechanism (the spur model^^) became very widespread in spite of ignorance of 
the presence of positive molecular ions and other intrablob electrons, different mobilities of e+ 
and e~ (Sec. II). Anyhow, combination of the Ore model and the spur model satisfactorily 
reproduces experimentally observed variations of Pps vs density and electric field. 

Since the mid-1970's the Ps formation at low fields (< 200 kV/cm) was also investigated 
in liquid and solid paraffins^^'^° and some pure liquids. ^^^^^ Below several tens of kV/cm the 
Ps yield noticeably drops, then the falling decelerates and a tendency to reach a plateau is 
observed. Application of much stronger fields (~ 1 MV/cm) may push the positron in the Ore 
gap and keep it there until the Ps will be formed, or e"*" annihilate with one of the molecular 
electrons. Thereby the strong electric field stimulates reaction (1). It is very probable that this 
phenomenon was experimentally observed in solid (polar and nonpolar) polymers^^, though 
the other factors (formation of bulk and surface charges, field induced electron and positron 
trapping/detrapping on structural defects) were mentioned, which may also increase the Ps 
yield in the fields about 1 MV/cm.^^ 

A theoretical interpretation of these experiments in condensed phase was attempted by 
Brandt et al.^^'^'' basing on the theory of Teutsch and Hughes^ developed originally for gases. 
The theory successfully describes the growth of Ps yield in strong fields because of the "heating" 
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of the positron and keeping it in the Ore gap. Characteristic field strength (when the Ps yield 
begins to grow up) may be estimated from the condition LET = ei?*, where LET (the linear 
energy transfer of e"*", —dW/dx) is determined by excitations of molecular vibrations (typical 
value is about 10"^ eV/A). Thus we find that the corresponding field is £J* ~ 1 MV/cm. 

In small fields the situation seems to be different. According to the Brandt theory, the 
decrement of the Ps yield is explained by the influence of the fleld on the positron escaping 
from the Ore gap during its thermalization. This process is considered as e^ random walks in 
energy space with the field dependent diffusion coefficient Dp{E) ~ {eElY /{l/v). It is expected 
that, when the field is applied, the positrons should abandon the Ore gap faster and therefore 
Ps formation must decrease. However, it is unclear why the authors considered only positron 
diffusion out of the Ore gap but did not include their diffusion back to this energy region. 
This theory completely ignores the track effects on Ps formation. It is unable to interpret 
the suppression of the Ps yield to almost zero when a quasifree electron scavenger is added 
to the medium. These effects are beyond the Ore model and the theory^^'^'''. Moreover, it 
predicts too large thermahzation lengths for e+ and e~, 100 times higher than ones known in 
radiation chemistry.^® As we shall see in the next sections these effects in low fields get a natural 
explanation within the framework of the recombination mechanism of the Ps formation. 

In the present paper we have developed new approach for consideration of the electric field 
effect on Ps formation at low electric fields (< 100 kV/cm), where it is possible to neglect 
positron acceleration by the field up to the eV energies. This approach is based on the recom- 
bination mechanism (diffusion-recombination model or blob model), ^^^^^ and properly takes 
into account the multipair nature of the end part of the positron track. 

In Section II Ps formation in condensed molecular media is considered. We introduce a 
concept of the quasifree positronium (a weakly bound e"'"-e'-pair, a precursor of the Ps in 
a bubble). We also present the arguments in favor of that the Ore process plays negligible 
role at low electric fields, where an adequate description of Ps formation may be achieved in 
the framework of the recombination mechanism. Finally the difference between spur and blob 
models is discussed. In Section III an approximate mathematical formulation of the problem 
is given. General expression for the Ps formation probability is obtained and some particular 
cases are considered. The last section contains some discussions of the results and further 
extension of this model. 

II. POSITRONIUM FORMATION IN CONDENSED MEDIA 

A. Quasifree Ps state. Modification of the Ore gap 

Consideration of the Ps formation in condensed molecular media (dielectric liquids, poly- 
mers, some molecular crystals) requires clarification of all intermediate stages and states, pre- 
ceding the formation of the final equilibrium state of the Ps atom. It is well known that in 
liquids the repulsive interaction resulting from the exchange of the electrons between Ps and 
surrounding molecules eventually leads to the formation of an equilibrium bubble with the 
radius about some angstroms. However, as we shall see below, formation of the quasifree Ps 
(qf-Ps) precedes the formation of the Ps bubble. The qf-Ps state corresponds to the bottom 
of the lower energy band available to the interacting e"'"-e~ pair before any rearrangement of 
molecules takes place. The notion of the qf-Ps was invoked in the frameworks of the recombi- 
nation mechanism of the Ps formation in order to explain the Enfield effect on Ps formation in 
hydrocarbons^^ and in relation to the Ps bubble model. 
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An availability of a large free space in gases always allows neglecting the zero-point kinetic 
energy of Ps, caused by the presence of the gas molecules. So in this case Ps binding energy 
is simply —Ry/2. In condensed phase the presence of molecules, firstly, essentially increases 
Ps zero-point kinetic energy arising owing to Ps repulsion from them (both e~ and e"*" are 
repelled from the cores of atoms because of exchange and Coulombic repulsions, respectively). 
Sometimes it is called as a "confinement" of the Ps. Secondly, molecular electrons screen 
the e'^-e^ Coulombic attraction and increase average e"''-e" separation in the positronium. 
Obviously these factors decrease binding energy between the positron and electron constituing 
the positronium and therefore reduce the width of the Ore gap, which may even completely 
disappear. 

For the first time the confinement of the Ps was studied by Brandt"''^ on the base of the 
solution found by Sommerfeld and Welker^^'"^^ of the quantum-mechanical problem of a hydro- 
gen atom confined by an impenetrable sphere of the radius R. Strict analytical solution of the 
same problem for the Ps atom is not available because of the impossibility of formulating the 
problem in terms of the center-of-mass coordinate and relative coordinate when the external 
potential is present. 

Elementary estimation of the Ps confinement effect is the following. Ground state energy of a 
point particle with the mass 2m in an infinite potential well of the radius R is = ^ (^^) > 
where as = 0.53 A is the Bohr radius. For R = -nas ~ 1.66 A this energy becomes equal to 
the binding energy, —Ry/2, of the free Ps atom taken with an opposite sign (for H atom it 
happens at i? = 1.835aB = 0.97 A). This gives a hint, that Ps might become unstable with 
respect to a break up on e+ and e~ and indicates that in a dense medium the width of the Ore 
gap becomes narrower and even may disappear. 

Certainly such an approach simulating Ps repulsion from molecules by an infinitely deep 
potential wall is crude. A more realistic evaluation of the energy of qf-Ps in a liquid phase can 
be obtained basing on the modified Ps bubble model, ^^'^^ particularly using the values of the 
depth U of potential wells of equilibrium Ps bubbles {U does not depend on the radius of the 
Ps bubble). 

The matter is that the physical meaning of U coincides with that of the Ps work function 
V^^*^ (the energy needed for Ps to enter the liquid without any rearrangement of molecules and 
stay there in the quasifree state). With the free- volume radius of the bubble tending to zero, the 
positronium transfers from the bubble just to the qf-Ps state, which has no preferential location 
m a bulk. Note, that V^"" is an analog of the work function Vq of the quasifree (excess) electron, 
Cqf, but there is one important distinction. Vq~ is at the same time the ground state energy of 
the e^f, because the energy of the electron at rest after having been removed from the liquid 
to infinity is defined to be zero.^^ Ps work function Vq^ differs from the qf-Ps ground state 
energy, V^^^ — Ry/2, by a constant shift, because its energy far outside the liquid is not zero, 
but —Ry/2. However, it is convenient definition because, would we consider dissociation of Ps 
in vacuum, the energies of the largely separated e+ and e^ become equal to zero. 

Energetics of the qf-Ps state may be illustrated with the help of the Born-Haber cycle: 
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Gas phase: e+ + e" — - — > Ps 



/ 


k 1 

-V 


-Vo" 








p ^ 





Liquid pliase: e+„f + e" f qf-Ps 



It consists of the following blocks: 1) removal of the quasifree electron and positron from the 
liquid to vacuum accompanying with some expenditure (or release) of energy, namely the sum 
of the respective work functions, —Vq~ — V^^, taken with negative sign; 2) formation of the Ps 
with the binding energy —Ry/2 in vacuum; 3) sending back of the Ps atom into the liquid 
phase supplying it with the energy Vq^. From the Born-Haber cycle we are able to estimate 
binding energy, Ef,, of the qf-Ps state: 

Eb = -V,- - + C - %/2 < 0. (4a) 

It is believed that V^^ is negative because of prevalence of the polarization interaction of e^j 
with the medium molecules over repulsion from the nuclei. So we may roughly equate Vf^ to 
the energy P„ of the polarization interaction of e~£, which can be calculated by decomposition 
of Vq~ into the sum of kinetic K and polarization P_ terms: = K + P^.^^ For hydrocarbons 
experimental values of Vq" are rather close to zero, which means that K and P_ tend to cancel 
each other. Using relevant expressions from^® we obtain that V^^ ^ ^ —{2 to 3) eV. 
Further, from the experimental data on the ortho-Ps lifetimes and the widths of the "narrow" 
component of the 27-angular correlation spectrum accounting curvature dependence of surface 
tension, it is possible to obtain U = V^^. It was found, that practically in all investigated 
liquids of various chemical nature values of V^^^ are close to 3 eV. Thus from Eq.(4a) we obtain 
that Ef, ^ —(0.5 to 1) eV. It means that qf-Ps is a loosely bound structure although the 
long-range Coulombic attraction between e"*" and e~ always provides an existence of the bound 
state (at least within the Onsager radius). Because positron-electron separation, Vep, in qf-Ps 
is expected to be large we may rewrite Eh in the following form: 



(46) 

2 ei 



where Eoo is the high-frequency dielectric permittivity of the liquid, and 3aB is the distance 
between e"*" and e~ in the Ps atom in vacuum (Ps diameter is {ip\rep\(p) / {ip\ip) = Sqb, where 
(p oc exp{—^rep/aB) and /i = 1/2 accounts for the reduced mass factor). The scaling E^ oc 
directly follows from the Schrodinger equation for Ps and is confirmed experimentally.^^ Using 
above mentioned numerical assessment for Ei, and Eq.(4b) we conclude that r^p is about 5 A. 
In such a Ps state delocalized over r^p the density of the "native" electron on the positron is 
rather small and e^ will primarily annihilate with electrons of the medium with the average 
lifetime about 0.5 ns and contribute to the free e"*" component of the AGAR spectrum. Thus, qf- 
Ps component could hardly be resolved experimentally, but it should renormalize contribution 
ascribed to the "free" e"*" annihilation. If the standard three exponential decomposition of the 
LT spectra is used the presence of qf-Ps is revealed, for example, in anomalous behavior of free 
e"*" lifetime.^° 

Further evidence that qf-Ps is a delocalized object is provided by rather high values of 
y-Ps pi2;g|;iy^ jf accept for the moment that qf-Ps in liquid may be considered as a "point" 
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particle, in that case Vq should represent just its zero-point kinetic energy due to the repulsion 
from surrounding molecules. Secondly, it is quite reasonable to assume that an excess quasifree 
electron injected into the liquid undergoes the same exchange repulsion as the " point" -qf-Ps. 
Therefore, it immediately implies that zero-point kinetic energy, K, of the e~£ should be twice 
as large as that of positronium and equal to 2V^^ ~ 5-7 eV. It is simply because the electron 
is half as heavy as Ps. However, it is very difficult to conceive that the negative polarization 
contribution, P_, could compensate 2Vq^ term and yield known experimental values of Vq~ , 
which are usually close to zero or even slightly negative. 

Finally we are able to estimate the low boundary Wiow of the Ore gap in condensed phase. 
Because the Ore process is just an electron transfer reaction, we assume that no rearrangement 
of molecules occurs and, therefore, final positronium will be in the quasifree state (formation 
of the bubble requires much longer time). If translational kinetic energy of the final qf-Ps 
is thermal (~ 0), the excess kinetic energy, W, of the projectile positron in this case should 
coincide with the low boundary the Ore gap: W = Wiow Now, let us carry out the Born-Haber 
cycle again, i.e. virtually pull out reagents to the gas phase, accomplish Ps formation reaction 
(1) and return products to the liquid: 



/e-Ry/2 

Gas phase: e+* (W) + M > Ps + M+- 
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In parenthesis we displayed the total energy of the e"*" in the liquid and gas phases. AG^ap is 
the vaporization energy of the molecule (a few tenth of eV) and P+ (negative) is the energy of 
polarization interaction of the radical-cation M"*"' with environment. Calculating total energy 
balance, one obtains that Ps formation may proceed if + — AG^ap — {Ig — Ryf^) ~ 
Vq^ + |P+| > 0. Using an approximate relationship between the liquid phase and the gas phase 
ionization potentials, II = Iq — \P+\ + Vg", it is possible to write the following relationship for 

Wio^ = h- V - V + ^G^ap + V^' - Ry/2. (5) 

Adopting that — V^q" ^ 9 eV (Table I) and taking into account previous estimations for other 
quantities involved, we obtain that Wiow is 7-8 eV. It is close enough (or may be even higher) to 
the threshold Wex of electronic excitations. Thus, in molecular media the Ore gap Wex — Wiow 
is rather narrow or may completely disappear, Wex < Wiow- Smallness of Wex — Wiow together 
with rather high value of the e+ energy losses (~ 10~^ eV/A for subionizing particle)^^'^*^ related 
mainly to excitation of molecular vibrations, make e"*" residence time inside the Ore gap very 
short (less and probably much less than 0.1 ps). That is why, generally speaking, Ps formation 
according to the Ore mechanism is ineffective in molecular media. However, it appears to be 
important when high external electric fields are applied. 
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B. Recombination mechanism of Ps formation. Spur and blob models 

Recombination mechanism postulates that Ps formation proceeds via combination of the 
positron and one of the knocked out electrons in the terminal part of the positron track. There 
arc two models, which utilize this mechanism, namely, the spur model^'^^ and the blob model 
(or diffusion-recombination modcl)^'^^'^^'^^. According to these models the decrease of the Ps 
yield in low fields (up to 100 kV/cm) occurs because the field takes apart the positron and 
track electrons in different directions. Therefore, probability of their encounter, as well as the 
Ps yield, drop down. In spite of the fact that both models answer the question about the Ps 
precursor in the same way, they differ in the point what is the terminal part of the e"*" track 
and how to calculate the probability of the Ps formation. Quantitative formulation of the spur 
model was given by Tao.*^^ It is based on the following assumptions: 

1) the positron and secondary electrons knocked out during of the last several acts of 
ionizations thermalize within the same volume of the medium. Positive ions also reside in the 
same space; 

2) secondary electrons recombine with the same probability cither with the positive ion or 
the positron. Therefore, the probability that the positron gets an electron for subsequent Ps 
formation is equal to tiq/ (1 + no), where Uq is the "initial" number of ion-electron pairs in the 
terminal spur; 

3) probability of the "elementary act" of e"'"-e~ recombination is taken in the Onsager form, 
l—e'^"/^^ where Tc = e^/sT is the Onsager radius, T is the temperature in energy units, b is the 
"initial" distance between e"*" and e~ by the end of their thermalization, and e is the dielectric 
permittivity of the medium. 

Thus, in the frameworks of the spur model Ps formation probability is written as: 

Pps = - e-'-^/^) • e-^^*p=. (6) 

1 -|- 77-0 

The last exponential factor takes into account a possibility of the free positron annihilation 
during the Ps formation time tpg (some picoseconds) with the annihilation rate A/ (< 2 ns~^). 
Obviously, contribution of this factor is negligible. 

Because in all nonpolar molecular media e ^ 2 at room temperature Tc is equal to ^ 300 
A. On the other hand, radiation chemistry data''^ show that typical thermalization lengths 
b of electrons are < 100 A. Thus, the Onsager factor, 1 — e~^=/^, is also very close to unity. 
Therefore, to explain observable values of Ps yields, which never exceed 0.7, we have to adopt 
in accordance with Eq.(6) that the terminal positron spur contains in average 2-3 ion-electron 
pairs. 

Later on, the recombination mechanism of Ps formation became extremely widespread. It 
is able to interpret extensive data on Ps chemistry (see^ and ref. therein). In contrast to the 
Ore model, this mechanism easily explains variation of the Ps yields from up to 0.7 in the 
substances of various chemical nature, where the Ore gap is practically the same. Changes of 
the Ps formation probability under phase transitions also received natural explanation. Exper- 
imentally observable monotonic decrement of the inhibition curves of the Ps yields (practically 
down to zero) in solutions of electron acceptors contradicts the Ore model, but is well incorpo- 
rated to the recombination mechanism. It naturally explains antiinhibition and antirecombi- 
nation effects, including experiments on formation Ps in low electric fields in pure liquids and 
mixtures. ^^''^^'■^ 

In*^^ the spur model was used as the basis for computer simulation of the Ps formation 
process. Terminal part of the e"*" track was approximated by a linear sequence of separated 
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ion-electron pairs. Fitting of the data in 2,2-dimethylbutane, it was obtained that the average 
distance between final acts of ionizations at the end of the e"^ track is about 165 A, which 
seemingly agrees with the spur model. Nevertheless, having considered the data for more 
than 50 liquids, the authors^^ came to the conclusion that multiparticle effects in the terminal 
positron spur are rather important for correct description of the processes in it. 

However, inconsistencies of the spur model should be mentioned (and the conclusions of the 
above simulations, in particular) with the data on the structure of the positron track (Fig.2). 
The matter is that the end of the positron track is not a spur, containing 2-3 ion-electron 
pairs, but the blob, containing about 30 overlapped ion-electron pairs. In essence this 
statement relies on the following two well known facts related to behavior of the energized 
electron (positron) on the blob formation stage, i.e. at energies Ry '^W < 500 eV (Appendix 
A): 

1) Ionization energy losses of such a positron (electron) is about 1-4 eV/A. It means, that 
e"*" very effectively ionizes molecules at the end of its track. Because the average energy Wiep, 
needed for creation of one ion-electron pair is^ i/Giep ~ 16-22 eV, (Gjep is the initial yield 
of ion-electron pairs per 100 eV of the absorbed energy) the distance k between subsequent 
ionizations made by e"*"* is about Wiep/hET fti 10-20 A. Dependence of k vs e"*" energy is shown 
in Fig. 3; 

2) Motion of the positron in the blob is diffusion-like (its transport path at these energies 
varies from 20 down to several A, Fig. 3). Thus, efficient e"*" scattering at the blob formation 
stage results in approximately spherical shape of the blob. Therefore, the simulation of the 
end part of the e+ track by the straight-line sequence of isolated ion-electron pairs'^^ or by a 
cylindrical column of ionizations^^ seem to be unjustified. 

The statement that the terminal part of the e"*" track is a blob, but not a spur is not just a 
question of terminology. As we shall see below, some primary processes, including Ps formation, 
proceed in a different way in blobs and spurs. Another important distinction is that in the blob 
model behavior of the positron is considered as very different from that of intrablob electrons 
and ions. Usually in the blob approach it is assumed that e+ is rather mobile and easily escapes 
from its blob (Sec.II E). Particularly, it imphes that the multiparticle nature of the terminal 
part of the e"^ track cannot be correctly taken into account via the factor no/ (no + 1)- In details 
we shall discuss these differences below, but now we shall briefly remind the life history of the 
positron, starting from its birth up to the annihilation. 

C. Life history of the positron 

Even though the positron is a stable elementary particle, its lifetime in matter turns out to 
be very short owing to presence its antiparticles - electrons. Usually the positrons are born in 
nuclear j3~^-dccsij (of ^^Na or ^^Cu) with typical initial energy about several hundreds of kcV. 
After implantation into a medium their life may be divided into three stages. It is convenient to 
subdivide the first stage into two steps: a) slowing down by ionization and b) thermalization. 
Retardation from the initial energy down to the ionization threshold lasts about 10~^^ s. When 
e+ energy becomes less than ~ 500 eV, the positron forms its terminal blob (Appendix A) . 
In the second step, the subionizing positron and knocked-out electrons become thermalized, 
exciting primarily intra- and intermolecular vibrations. In molecular liquids, it usually takes 
less than 10~^^ s in addition to the duration of the first step. 

The second stage involves fast intrablob reactions with the participation of primary radical- 
cations RH+" (ions) and thermalized electrons and positrons. According to the recombination 
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mechanism, Ps is formed just at this stage. This stage is hmited by locahzation (solvation), 
ion-electron recombination, or out-diffusion of the species. The typical duration of this stage 
is 10-^2_io-ii s. 

The third stage includes different chemical reactions between Ps, locahzed species, and 
additives. It is terminated by e"*" annihilation from para-Ps, "free"-e'^, or ortho-Ps by the 
pick-off process. This final stage lasts about 10~^°-10~^ s. 

D. Main distinctions between spurs and blobs 

1. Sizes and energy release 

In radiation chemistry the term " spur" denotes the small isolated spherical part of the track 
of an ionizing particle, containing a few (tIq^ = 1-5) ion-clcctron pairs. ^^'^^ Average energy 
release in the spur is 30-50 eV. Its typical size is determined by a thermalization length of the 
knocked out electron in the Coulombic field of the ion. In a crude description of radiation- 
chemical processes it is possible to assume, that tracks of high energy electrons and Compton 
electrons produced by ^'^Co 7-rays consist only of isolated spurs (contribution of short tracks 
and blobs is rather small). In spurs with a few ion-electron pairs, Onsager's description of the 
geminate recombination is plausible. 

As we already mentioned, the end part of the e"*" track is a blob. Calculations show (Ap- 
pendix A) that the spatial distribution of ionization events in it can be simulated by gaussian 

function no— ^-^sttt"^- Here Uq — W^/Wiep ~ 30 is the initial number of ion-electron pairs and 

^ bl 

40 A. The blob size is determined by the ionization slowing down of the energetic 
positron, when it loses energy from ~ 500 eV (total energy deposition in the blob) down 
to the ionization threshold. 

2. Quasi-neutrality condition 

In blobs knocked out electrons before getting localized rapidly adjust their motion to the 
distribution of the primary positive ions and efficiently screen them^^. Such a quasi-neutrality 
electron distribution in the field of ions is established because 1) Coulombic interaction between 
ions and electrons is strong enough (the Onsagcr radius is much larger than the blob size and 
average electron thermalization displacement"^^) 2) mobility of electrons is much higher than 
that of ions. Simple estimations (Appendix B) show that the width of the spatial distribution of 
the intrablob electrons is only a few angstroms larger than that of ions: ag — Oj ~ 2.4a^;/rcno -C 
Qbi. Hence ^ ^ Ui. Here and af are dispersions of the distributions of the electrons 
and ions, respectively. The concept of quasi-neutrality of the blob seems reasonable especially 
in non-polar liquids where an absence of significant spatial fiuctuations of local electric field (in 
contrast to polar media) favors establishing blob quasi-neutrality and complete recombination 
of charges. 

3. Expansion of the blob 

Because of the Coulombic attraction between charged intrablob species their diffusion expan- 
sion does not proceed independently. As a result out-diffusion of electrons is almost completely 
suppressed, but the diffusion coefficient of ions is increased by a factor of 2. This case is known 



9 



as ambipolar diffusion (Appendix C). Thus blob expansion proceeds very slowly and may be 
practically neglected in the problem of Ps formation. 

4- Elongated thermalization and recombination 

After the last ionization event in the terminal blob, the positron and intrablob electrons 
become subionizing. Having no possibility to excite electronic transitions in molecules, they 
lose energy by exciting molecular vibrations. In polar media there is an additional way of 
energy losses: reorientation of the dipole molecules (Debye losses). Thermalization theories 
usually consider molecular excitation a result of its perturbation by an alternating electric field 
induced by a flight of an isolated charged particle. ^^'^^ If the Fourier spectrum of this field 
contains harmonics with frequencies close to those of molecular vibrations, excitations get very 
probable. In the blob because of the presence of a large number of ion-electron pairs and 
efficient screening (local quasi- neutrality) , the resulting electric field at a given point does not 
fiuctuate significantly. Thus, in spite of fast motion of intrablob electrons, they cannot lose 
kinetic energy easily via excitation of vibrations. This lengthens thermalization of electrons and 
ion-electron recombination because these processes depend on the absorption rate of released 
energy. 

E. Behavior of the positron and positronium formation 

1. e"*" escape from the blob 

Because the blob is electrically neutral, Coulombic interaction between the positron and 
its blob is absent. Polarization attraction between c'*' and the blob is small enough (Sec.V 
F). Therefore, after the last ionization even subionizing positron may easily to escape from its 
blob during thermalization by travelling a longer distance than the blob electrons, which are 
held by the ions. So one may expect that the width of the e"*" spatial distribution by the end 
of thermalization is larger than the blob size: Op > ai,i. Strictly speaking is a dispersion 
of initial (gaussian) distribution of the thermalized positron around the center of the blob. 
Because e"*" thermalizes primarily outside the blob its thermalization time is not affected by the 
above mentioned elongation of the thermalization process. 

2. e"*" mobility 

In non-polar liquids the positron after its thermalization is expected to be in the quasifree 
state having rather high mobility in comparison with than that of the quasifree electron (Table 
I). Firstly, it is due to the absence significant spatial fiuctuations of local electric field, which 
may trap or localize the positron (the medium is nonpolar). Secondly, a density of a liquid 
state is too high to permit clasterization of molecules around the positron like in pressurized 
gases. 

Actually because of the screening, Coulombic interaction with nuclei plays minor role in 

the scattering of the quasifree positron and electron. Their mobilities are primarily determined 
by the small angular scattering on outer molecular electrons (see, for example, Eq.l2 in ref.^°). 
However, for excess electrons there is some spin-exchange contribution which enhances e~ scat- 
tering and decreases their mobility in comparison of the ei one. Exchange repulsion may even 
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lead to full localization of the excess electron in a void. Contrary, the positron prefers to reside 
in a bulk because of polarization interaction. Preferential e^ traps are positive density fluctua- 
tions (aggregation of molecules). However, such traps practically cannot be formed in a liquid 
phase because of strong asymmetry of intermolecular interaction potential. 

3. "Elementary act" of Ps formation 

Being thermalized primarily outside its terminal blob, quasifree positron may diffuse to 
the blob again and pick up one of the intrablob electrons. That process needs the "initial" 
separation between the particles not exceeding the average distance between intrablob species, 
i.e. (47ra^;/3no)^/^ ~ 20 A. According to Eq.(4b) initial binding energy of this e+ • • • e~ pair is 
about 0.1 eV. Translational kinetic energies of each particle composing the pair must be less 
than the initial binding energy, otherwise the pair will break up. Therefore e+ and e~ just 
before the "elementary act" of Ps formation must be practically thermalized. 

In the framework of the adiabatic approximation (the positions of molecules are fixed) such 
a pair is not at the bottom of its energy spectrum and may continue to release energy via 
excitation of intra- and intermolecular vibrations and finally form at equilibrium (in adiabatic 
sense) the qf-Ps state discussed in Sec. II. A. In comparison with the binding energy of qf-Ps, 
the e"*" ■ ■ ■ e~ pair possesses about 1 eV of excess energy, which is accumulated in the form of a 
potential energy of the mutual Coulombic attraction. Let us now see how this pair transforms 
into qf-Ps and what happens with this excess energy during e"^ • • • e" ^ qf-Ps transformation. 

Let total momentum of the e"*" • • • e~ pair be zero at the beginning. When particles getting 
closer, kinetic energy of both the positron and electron increases. If then one of the particles 
creates a phonon (excites a molecular vibration) its energy drops down (as well as momentum 
and it changes direction also) and as a result the total momentum of the e"*" • • • e~ pair becomes 
non zero (just due to momentum conservation in inelastic phonon scattering). This mechanism 
gives rise to the motion of the e"*" • • • e~ pair. Adopting ^ 0.01 eV/A as a typical value of the 
e~ (e"*") LET at energies about 1 eV due to the excitation of vibrations, we may estimate total 
path of the e"*" ■ ■ ■ e~ pair before it becomes qf-Ps. It is E'^/LET ^ 50 A (the value of Ei, is 
estimated in Eq.(4a)). When the positron and electron get closer, their polarization interaction 
with the medium decreases. Larger number of surrounding molecules become to "feel" e"^ • • • 
e~ pair as an electrically neutral object. So LET of such a pair may decrease and its total path 
may increase. 

After losing all its excess potential energy the e"*" ■ ■ ■ e~ pair becomes the qf-Ps, which 
is in thermal equilibrium with the environment. In a disordered medium (for example, in a 
liquid) it immediately finds nearest place with slightly lower density and will stay there, if the 
temperature cannot push it out from that shallow trap. This state may survive for a rather 
long time (it may be comparable with the free e"^ lifetime) in a solid medium, where molecular 
rearrangements are energetically forbidden or require longer time. 

In liquids (or "soft" media) qf-Ps state continuously evolves into the bubble state at a 
longer time scale because of further gain of the e^-e" Coulombic attraction and decrease of the 
positronium zero-point kinetic energy. It becomes possible due to rearrangement of molecules 
and appearance of additional free space. This is the beginning of the Ps bubble formation stage. 
Kinetic energy released during compactification of the Ps stimulates further local decrease of 
the density. Eventually an equilibrium bubble with the Ps inside is formed. Schematically Ps 
formation may be displayed as follows: 

ejf -I- e^Qb — > e"*" • • • e" — > qf-Ps Ps in the bubble. (7) 
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Cursory consideration of the bubble growth based on the solution of the Euler equation 
for incompressible fluid taking into account the Laplace pressure and the quantum-mechanical 
"pressure" exerted by the Ps atom on the wall of the bubble, gives the following estimation for 
the bubble formation time 



where p is the density and a is the surface tension coefficient of the liquid. It implies that 
bubble formation lasts about several picosecond in the majority of liquids. On the similar basis 
consideration of the self-locahzation of the positronium (the last step of Eq.(7)) was suggested 
by Khrapak.^^ 

Perhaps molecular crystals (anthracene, acenaphthylene) give us the best example illustrat- 
ing the transformation of qf-Ps into the bubble state. It is known that in these substances in 
the solid state there are neither longlived nor "narrow" components in the lifetime and angular 
correlation (AGAR) spectra, respectively. So from a viewpoint of experimentalists the Ps atom 
does not form there. However when the temperature is slightly above the melting point (in the 
liquid phase) Ps formation probability becomes very large (more than 0. 6 '"''^ According to 
our model, it is due to impedance in the solid phase of the qf-Ps transformation to the bubble 
state, but in the liquid phase it easily proceeds. As a result, in the solid phase we expect con- 
tributions both from free e"*" and qf-Ps annihilations, but it might be difficult to resolve them 
experimentally. Annihilation parameters of qf-Ps (lifetime and width of its AGAR component) 
may be obtained by subtracting from the solid-state annihilation spectra the free e"^ contri- 
bution, which must be taken (with the same intensity and lifetime; corrections accounting for 
the variation in the average electron density at the phase transition may be needed) from the 
measurements in the liquid phase. 



From a viewpoint of the energy balance localized electron and positron might participate 
in Ps formation. Actually Cj"^ already resides in a small void and its binding energy there is 
about 1 eV. Hence, Ps formation may proceed through e^^ trapping by e^^^ on an appropriate 
Rydberg-type orbit with a subsequent formation of the equilibrium Ps bubble. 

However, in non-polar liquids, as we shall see below, electron localization takes much longer 
than Ps formation and, therefore, ej"^ does not exist on a time scale of Ps formation process. 

In polar liquids, like water, one may expect that e+ gets also localized within a time compa- 
rable with electron localization. After that, mobility of the particles drastically drops down and 
they simply have not enough time to meet each other (time is limited by a free positron lifetime 
~ 0.5 ns). Really the diffusion displacement of ej|^^ during 0.5 ns is smaller than thermalization 
displacement (~ Cp). 

Thus, in what follows we do not include localized particles into equations describing Ps 
formation. 

III. EARLY INTRABLOB PROCESSES: GENERAL FORMULATION 

Thus, according to the blob model we assume that Ps formation proceeds via recombination 
of the thermalized but not yet localized positron with one of the intrablob quasifree electrons, 
released while the positron is being slowed down by ionization: 




4- Role of localized species 
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+ e~ — > qf-Ps — >■ Ps in the bubble. (8a) 
This reaction may compete with other intratrack processes 

locahzation : e e^^^, e^^ ej^^, (86) 

ion-electron recombination (lER): e~ + RH"'"' RH* — > Products. (8c) 

ion-molecule reaction: RH"*"' -|- RH ^ R' -|- RH^, {8d) 

In CgFg intratrack electrons by the end of their thermalization can be weakly captured by 
host molecules:^^'^^'^^ 

e- + C6F6^C6F6. (9a) 
However, the positron may easily pick up such an electron from CeFg ion and form Ps: 

<f + CeFg- ^ Ps + CeFe. {9b) 

Therefore in this case Ce(r, t) in Eqs.(lO) will stand not for the concentration of the electrons 
(as in other liquids), but for the concentration of CeFg anions. 

The blob model describes nonhomogeneous kinetics of reactions (8) in terms of concen- 
trations of the particles, because the number of particles involved is large and local motion of 
quasifree electrons and the positron is fast enough in comparison with Ps formation time (Sect.V 
E). In the simplest case the external electric field enters diffusion-recombination equations (10) 
via a drift of the positron distribution (electrons are tightly bounded to ions) : 

dc^{r, t) ^ _ kieCiCe - Ci/Ti, (10a) 



dt 

dt 
dcp{v,t) 



— ^amb^^e kigC-iC^ k^pC^Cp C^/t^, (106) 

= DpAcp - div(t;pCj,) - kepCeCp - Cp/Tp. (10c) 

Here Ci{r,t), Ce(r, t), Cp{r,t) are the concentrations of the primary radical-cations RH+' (ions), 
intrablob electrons and the positron, respectively, at time t at point r (measured from the 
center of the blob), A is the radial part of the Laplace operator, Damb = '^Di is the coefficient 
of ambipolar diffusion (Appendix C), Di and Dp denote the diffusion coefficients of ions (i) and 
the positron (p), respectively, and hie and k^p are the recombination rate constants. Decay rate 
includes electron localization as well as possible capture by host molecules, accompanied 
by formation of their radical- anions (in this work we study pure liquids only). Similarly, is 
the rate of the ion- molecule reaction and l/r^ = l/r2 + ^I'T^p'^ accounts for the e"*" annihilation 
{t2 is the free positrons lifetime. Table I) and possible localization of the positrons with a 
characteristic time t^°'^. 

As discussed above and in Appendix B, blob species distribution is only slightly disturbed by 
the application of an external electric field, while the behavior of the positron is very sensitive 
to the presence of the field. The field decreases Ps formation by removal of the positron away 
from the blob, which suppresses reaction (8a). In Eqs.(lO) this effect is approximately taken 
into account through the term, div(i;pCp), describing the drift of the positron in external field 
D. However, we skip polarization interaction between the blob and the positron (see Sec. V G) 
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and neglect deviation of the field from its average value within the blob caused by the presence 
of ion-electron pairs. We also assume linear proportionality between the drift velocity Vp and 
D, i.e., Vp = bpFi = bpD/e {bp = eDp/T is the positron mobility and e is the static dielectric 
permittivity of the medium; see comments to Table I). This assumption is justified for D < 100 
kV/cm where the nonlinear effect is not significant. A more exact way to account for real field 
distribution at the terminal positron blob is a joint solution of Eqs.(lO) and an equation on the 
electric field E(r) as was done in Appendix C. 

Because the motion of the knocked-out electrons during blob formation and subsequent 
thermalization obeys the quasi-neutrality condition, it is natural to adopt the following initial 
conditions for Eqs.(lO): 

c.(M=0) = c.(M=0)=n,.^^P(3l!^. c,(M=0) = („) 

Hq is the initial number of ion-electron pairs and a^^ and are the dispersions of spatial 
distributions of the intrablob ion-electron pairs and the positron. These distributions do not 
take spatial electron-positron correlations of into account. In contrast with the blob electrons, 
the motion of subionizing positron is diffusion-like and it may easily escape from the blob. 
Hence, dispersion of the initial positron distribution may significantly exceed that of intrablob 
species. In contrast to ahi, Qp is determined by thermalization of the positron, i.e., by its ability 
to excite intra- and intermolecular vibrations, while the distribution of the ions in the blob 
(i.e. au) depends on a much more efficient energy loss process - slowing down of the energetic 
positron by ionization, when its energy is reduced from Wu to the ionization potential II- In 
Eq.(ll) Ip accounts for the drift of the subionizing positron during its thermalization: 

Ip = Vp{t)dt = ^ Dp{t)dt = ^ . (12) 

Here we used that the average e"*" thermalization displacement squared 6 jj*'' Dp{t)dt during e+ 
thermalization time tth is equal to — ali. 

To calculate Ps formation probability, Ppg, we must integrate the term kgpCeCp over the 
whole space and time: 

Pps = KpJ dt J Ce{r,t)cp{r,t)d^r. (13) 

This is the probability of formation of qf-Ps. However in liquids we expect that this state 
transforms to the bubble state rather fast, within several picoseconds as was discussed in the 
previous section. 

Within the framework of the prescribed diffusion method^^ the solutions Cj{r,t) can be 
written in the following form: 



{r,t)^njit)Gjir,t), J Gj{r,t)d'r ^ 1, j^{i,e,p}. (14) 



ne(t) and ni{t) are the numbers of quasifree electrons and ions survived up to time t, and np{t) 
is the free positron survival probability. Gj{r,t) is the Green function of a simple diffusion 
equation (without recombination and decay terms): 

_ exp[-(r - 1, - Vjt)y{4Vjt + a])] 



14 



We introduced different subscripts for Vj, Vj, aj and Ij simply to maintain symmetry. In 
accordance with Eqs.(lO) and Appendix C Vi ^ 0, Vg ^ 0, Vi = = Damb <^ T^p = Dp, 
tti ^ tte ^ au and Ij = Ig = 0. Substituting Eq.(15) and (14) into Eqs.(lO) and integrating the 
resulting equations over whole space, we obtain a much simpler system of ordinary differential 
equations on nj{t): 

rii = -kieneUi/Vie - rii/Ti, ni{0) = Uq, 



fie — —kieTleni/Vie — k^pTlgnp/Vep — Ue/Te, ^^-6(0) — Uq, 



rir. 



where 



GjGkd^r 



exp 



np{0) = 1, 

{\j + Vjt - Ifc - Vkt)'^ 

{a] + al){l + t/Tjk) 
j,k = {t,e,p}. 



(16a) 
(166) 
(16c) 

(17) 



In Eq.(166) for rie, it is possible to omit the term kepUenp/Vep, because it has a negligible 
effect on the disappearance of the intratrack electrons. Then within the prescribed diffusion 
approach Eq.(13) is rewritten as follows 



~ ^ dt — kf.p I 

Vep Jo 



^e^p 



V^p{l + t/repr/' 



exp 



2 _ ya + t/Tep) 
ep ' 

' ep 



1 + t/r, 



dt, 



(18) 



where 



eD 



(^bl + 



I ep 



' ep 
oil + «r 



AeT 



2 d. 



3 al + al ■ 



If we neglect positron annihilation and localization, i.e. put ^ = 0, one obtains 



Pi 



Ps 



(19) 



np(oo) is the fraction of Ps formation and 



This equation is a conservation law for positrons: 1 
the remaining np{oo) is the fraction of positrons annihilating in a "free" state at a time scale 
on the order of T2 = Aj^. 

Application of the prescribed diffusion method to solve Eqs.(lO) is based on an assumption 
that interaction between thermalizcd positron and its blob is negligible. In other words, we 
adopted that the blob is "transparent" for e^j (wc call this approach the "white blob" model). 
This assumption seems reasonable because of two opposite effects which approximately cancel 
each other (Sec.IV.G): 1) outdiffusion of intrablob electrons makes them to reside in an outer 
region of the blob, which results in appearance of an excess positive charge in its central region, 
repelling e+ from the blob; 2) the presence of e"*" within the blob may lead to rearrangement 
of intrablob electrons, which may decrease the total energy of the system because of the Debye 
screening. 

Thus, in this case fast diffusion motion of the quasifree positron on a timescale of its lifetime 
(as we shall see in Sec.IV.E Tgp <^ T2) efficiently smears the positron distribution and approaches 
it to a Gaussian shape. This justifies application of Eqs.(14)-(15). 
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A. The simplest case of Ps formation. Relation to the Onsager approach 



If the Ps formation process is assumed to be very fast compared to lER, possible localization 
processes and diffusion expansion of the blob {kep ^ /cje — 0, Te — > oo, Tp — > oo and Dp ^ 
Damb): solutious of Eqs.(lO) become 



ne{t < Te 

where 



no, 



np{t < Tp) = exp 

^ok^pT^p 



-W, 



ep 



rt/Tei 

Jo 



(1 + 1?)3/2 



2 (a + 'i?)^ 



nok 



ep 



nok 



ep 



ep 



ep 



(20) 



(21) 



is a dimensionless parameter that integrally accounts for such factors as diffusion of the positron, 
e+-e~ electrostatic attraction, the efficiency of the absorption of released energy, free volume 
and its distribution, which determine the Ps formation probability in a zero field. Substituting 
Eq.(20) into Eq.(19), we obtain 



Pps = 1 - exp 



ep 



d-d 



{l+'&fl^ 



exp 



l-l'ep 



+ eD 



eT 



lO-\al + al{k) 



D /kV 



cm 



(22) 



(23) 



In the present experiments at room temperature, the highest experimentally reachable field 
was D = 25 kV/cm. If ^ a^^ + ~ 100 A, /Xep is about 0.1. It is worth mentioning that the 

drift VpTep — jJt'ep'sJ clIi + Op of the positron distribution during Ps formation time is small in 



comparison with the diffusion distance 



ttp during the same time. 



At a small fiep an asymptotic expression of Eq.(22) is: 



PpsilJ'ep < 1) = 1 - exp 



-2Wep { 1 - V^i^ep + ^(8 -5a-a^)-. 



(24) 



It is seen from Eq.(24) that a only slightly affects the Pps through the highest orders of /^ep- 
In a zero field (/Xep = 0): 



Pps(0) = l-exp(-2iyep). 



(25) 



Note that this equation reproduces the form of the well-known Onsager result for geminate 
recombination, which in low fields is written as:^° 



1-1 + 



eDvc 



exp - 



ro 



ir' 



(26) 



where tq is the initial separation of the geminate e"'"-e pair. If ^y/nWepfiep <^ 1, we can expand 
the exponent in Eq.(24) and keep only the leading field-dependent term: 



Ps 



eDvc 
2eT 



■ exp 



nok 



ep 



47rD. 



^^0 



n{ali + al] 



(27) 
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The meaning of radius in Eq.(27) differs from the Onsager radius used in Eq.(26). Eq.(26) 
corresponds to the 5-function distribution of e^-e~ initial separation, 5{r — ro), while Eq.(27) 
is obtained assuming Gaussian distributions for the positron and electron. 

Note that Eq.(26) contains only one fitting parameter, ro, while Eq.(27) has two. We thus 
have more freedom to fit experimental data with Eq.(27). With Eq.(26), we can fit only relative 
Ps yield Pps(-D)/Pps(0).'^'^''^'^^ This is an important difference between the Onsager model (spur 
model) and the blob model. It is also worth mentioning that the dielectric permittivity enters 
Eq.(26) and Eq.(27) in a different way. 

For a large /igp, the asymptotical expansion of Eq.(22) becomes 



PpsilJ'ep > 1) = 1 - exp 



ep 



1 - 



+ 



(28) 



2 Hep V 20r/i 

ep 

It means that Pps should go to zero at high fields. However, this regime is definitely beyond 
the applicability of the present theory (for example, inequality /lep > 5 implies that D should 
be larger than 10^ kV/cm). In this case another Ps formation process comes into play (Sect. 
Ill B). 



B. e"*" annihilation and ion-electron recombination 



A more realistic consideration than that in the previous section must take into account at 
least e"*" annihilation and the possibility of lER (reaction 8c). Both of these processes may 
compete with Ps formation. lER equally decreases the number of electrons and ions. If we 
neglect as before the diffusion expansion of the blob. 



ne{t) = ni{t) = 
In this case integration of Eq.(16c) gives 



no 



l + nokie{2nali)-yH' 



(29) 



np{t) = exp 







where 



'o (l + ^)3/2(l + W^,e^) 

_ npkieTep nokiejali + Op) 



2 (a + ^)^ 
exp ( -/iep 



(30) 



(31) 



Unfortunately, because of the presence of the decay term np/Tp in Eq.(16c) we cannot use 
the relationship Pps = 1 — np{oo) as before and must use Eq.(18) for calculation of Ppg. Its 
integration by parts gives 



Pps = 1 - Wp 



= 1 



(I'd exp 



ci'i?exp 



-'d-W, 



ep 



Wp=^-^ 



(^bl + (^p 



1 



1 

T2 T. 



loc 



(32) 
(33) 
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A distinctive feature of this relationship is a behavior at low fields. It is not a linear decrease 
with D as it was in Eq.(22), Eq.(24) and Eq.(27), but a quadratic one. At non-zero Wp or Wie-, 

straightforward expansion of the exponent exp ^— /^ep ^"i+z* ) ~ 1 — /^ep '^'^i+l small fields in 
Eq.(32) leads to 

Pp,(D = 0) - Pps(D 0) oc D^. (34) 

This important peculiarity is a result of competition between the Ps formation process, lER 
and e+ annihilation (or localization. 

IV. COMPARISON WITH EXPERIMENTAL DATA AND DISCUSSIONS 

A. Pps(£)): Simplest consideration 

Here we consider experimental data of Kobayashi and co-workers obtained by means of the 
positron annihilation lifetime spectroscopy in the following nonpolar liquids: benzene, hexaflu- 
orobenzene, hexane, cyclohexane, and isooctane.^^'^^ Air dissolved in the liquids was removed 
by the standard freeze-thaw method. High voltages were applied to liquids through a pair of 
electrodes to yield external fields of up to 25 kV/cm. In a given electric field, one lifetime 
spectrum was collected for 2.5 hours at room temperature (295 ±2 K), resulting in a total of 
~ 7 • 10^ counts for each spectrum. Measured lifetime spectra were decomposed into three 
components. The first and second components, ri =190-310 ps and T2 (Table I), were assigned 
to the annihilation of p-Ps and free positrons. The longest-hved component, T3, was due to the 
ortho-Ps (o-Ps) annihilation in bubbles in liquids (Table I). The o-Ps lifetime was unchanged 
by an increase in the external electric field. 

To compare our theoretical results with experimental data, we relate Pps to the intensity 
of o-Ps component I^. We assume that — 3Pps/4 • 100%. The multiplier 3/4 in this relation 
may vary slightly from one liquid to another because of possible Ps interactions with highly 
reactive radiolytic products (localized electrons, radicals, and radical-cations). ^'^'^ 

Experimental data arc shown in Fig.l. As a first step, we tried to fit experimental data 

to Eq.(22), adjusting two parameters: a^^ + and Wep (Fig.l and Table II). Parameter a 

(0 < a < 2/3) entering Eq.(22) was not free. It was recalculated from fi^p assuming = 40 
A. It is seen from Eq.(24) that a has small infiuence on Pps. In all liquids investigated but 

CqFq, ^ all + '^p 100-200 A. Obtained e~-e"'" recombination rate constants (Table II) are of 
order of lER rate constants measured in pulse radiolysis experiments in radiation chemistry: 
4.7 • 10^^-7.2 • 10^^ M-^s-^ in n-hexane, 1.9 • 10^^ M-^s"^ in cyclohexane and ~ 2 • 10^^ M-^s"^ 
in isooctane.^'' As mentioned, k^p inchidcs many different factors, e.g., c+-e~ electrostatic at- 
traction, availability of the appropriate free volume, and the absorption of the released energy, 
etc. 



B. Effect of the positron annihilation and ion-electron recombination 

It is seen that the electric field effect is stronger in saturated hydrocarbons (cyclohexane, 
hexane, isooctane) than in aromatic compounds. Experimental data in Fig.l in benzene and 

especially in hcxafiuorobenzene indicate that the slope of Pps(-D) tends to zero as the field 
approaches zero. This observation is reasonable if Ps formation competes with e+ annihilation, 
e"*" localization, or lER. Treatment of the data based on Eq.(32) (solid lines in Fig.l) with 
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Wie—0 and l/Tp°'^—0 leads to a zero slope of Pps(D) at D=0 and larger values of parameter 
all + '^p (Table II). The introduction of non-zero values of Wjg worsens the fit, but if we allow 

that Wie/Wep is 0.005, final curves coincide practically with solid lines, that is why we did 
not plot them. The values of fitting parameters are listed in Table 11. Based on these results, 
the maximal limiting value of the Wie/Wep ratio was found to be 0.005. Using the numerical 
values of 40 A and Uq ~ 30, the ratio of rate constants kie/kep turns out to be very 

small, about 10~^ (Table II), which implies that lER cannot compete significantly with Ps 
formation; note, however, that the values of kic estimated here do not correspond to typical 
experimental conditions in radiation chemistry. This indicates that, in the terminal positron 
blob, Ps formation proceeds within a shorter time than lER. The following three points make 
this clear: 

1. Both Ps formation and lER reactions depend strongly on the same process - energy 
absorption and the transfer of released energy. As discussed in Sect. II D.4, because of the 
high density of ion-electron pairs and quasi-neutrality condition in the blob, thermalization and 
lER are lengthened. Ps formation proceeds more easily, however, because both particles lose 
energy and the total released energy is less than that in lER. 

2. kie is proportional to the overlapping of wave functions of a delocalized state (quasifree 
electron) and a localized state (positive ion), while kep is proportional to a much larger over- 
lapping of the two delocalized wave functions of the quasifree positron and electron. 

3. As mentioned above the mobility of e^f is higher than that of e~f (Table I), leading to a 
larger ratio. Dp/ De, proportional to kep/ he- 



C. e+ localization 

To avoid the strong disagreement with experimental data, we must assume that e"*" local- 
ization does not proceed in neat liquids, t^"'^ at least should be longer than T2 (Table II), which 
agrees well with our previous conclusion (Sect. II E.2) that e+ localization is unfavorable in 
neat nonpolar liquids other than CqFq. 



D. CeFe 

The CeFe molecule may capture an epithermal electron in a shallow energy level. " Shal- 
low" implies that the positron may pick up an electron from CeFg and form a Ps atom. Trapped 
electrons escape lER but survive for Ps formation by remaining on a molecule of the liquid. 
Localization of intrablob electrons decreases the ■(/'-function overlapping with e^^, which results 

in a smaller rate constant kep (Table II). Product Dp^J ali + also turns out to be small, how- 
ever, Wep and therefore the Ps formation probability at zero field (Pps(O) ~ l — e^^^^'p, Eq.(25)) 
is the largest among the liquids investigated. The small value of the positron thermalization 
distance in CqFq qualitatively agrees with the data of Gee and Freeman for the thermalization 
distance of the e xcess ele ctron. 

Why is Dp ^ all + ^p small in hcxafiuorobenzene? The higher electron density Ue and the 

larger number Z of electrons in CgFg (Table I) enhance positron scattering and decrease Dp oc 
l/Zue- Thus e"*" mobility bp drops below the experimentally detectable limit. ^"^ 
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E. Ps formation time 



Using the numbers for y a^^ + and adopting approximate positron mobihty bp equal to 

10 cm^V^^s^^ (Table I), which corresponds to Dp = 0.25 cm^/s, we obtain Ps formation time 
Tgp ^ {ali + ap)/ADp on the order of some picoseconds for the hquids except CgFg in accordance 
with other estimations. ^^'^^ 

F. Absence of the field dependence of T3 

Experimental data clearly show that lifetime of the o-Ps atom residing in a bubble, does 
not depend on applied electric field within ranges of statistical deviations in all investigated 
liquids. It is well-known that the o-Ps annihilation rate is proportional to the overlapping of 
the positron wave function with that of the electrons of the molecules (" pick-ofF' process) . 

The Ps atom is electrically neutral, but a highly polarizablc system. Nevertheless an external 
field D ^ 25 kV/cm stretches e"*" and e~ in Ps only by a distance about ~ 10^^ A. It slightly 
increases the overlapping in one side of the Ps bubble and decrease it by the same amount on 
the opposite side. Total overlapping remains constant. This is the reason for the absence of 
the field dependence of T3 values. 

G. Electrostatic and polarization effects 

More careful inspection of our experimental data at high fields and other related experimen- 
tal results^"^'*^'^^ suggests that Ps formation probability tends to reach a plateau. It in fact could 
be attributed to a nontrivial role of polarization interaction between the thermalized quasifree 
positron and the blob, which is not yet well understood. 

Being highly mobile, intrablob electrons tend to reside in an outer region of the blob (Ap- 
pendix B). It results in the appearance of a small excess positive charge in its central region, 
which repels e"*" outward. This potential, Eq.(B7), is about 3T in the center of the blob and 
decreases to zero at r > a^. 

Polarizability au of the blob as a whole is about eafi, so at r > au, the e"''-blob polarization 
attraction is about — where r is e"*" separation from the center of the blob. Its maximum 

possible value (at r ~ au) is about — — ~ —0.1 eV, which is larger than T, so this shallow 
negative potential in the outer region of the blob may trap the thermalized positron. 

Because of the presence of the large number of ion-electron pairs inside the blob, e"*" potential 
energy decreases by a value ~ — — ytt- due to correlations in positions of the charged 

intrablob species, which is, in essence, Debye-Huckel screening of the positron charge. Here 
the Debye radius td ~ (47rrcCiep)"^^^ 4 A, where Qer ~ /fTra^^ is the concentration of 
ion-electron pairs in the blob, and au/nj^ is the average distance between intrablob species. 
At the distances less than the average distance between the particles r < au/ nj^ the screening 
potential is pure Coulombic. At larger r it takes the Dcbyc form. Such screening of the positron 
by blob charges also makes residence of the thermalized positron inside the blob favorable. 

These effects are rather subtle. Quantitatively, they are somewhat higher than T (~ 0.1 
eV), but have different signs in r-space. Nevertheless it is possible that e^f may be trapped 
inside or near the blob. This circumstance is important for interpretation of the electric field 
effect on Ps formation and may be related to a tendency of Pps reaching "plateau" at high 
fields. 
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The approach developed here properly takes into account the presence of the blob as an 
inhomogeneously distributed large number of ion-electron pairs and assumes that all of the 
above polarization effects more or less compensate each other. More accurate consideration 
of these electrostatic and polarization effects of the e+-blob interaction is the subject of our 
following study (so-called "black blob" model). 

V. CONCLUSION 

The difference between radiation chemistry and Ps chemistry is related to the difference in 
the objects they study. Being a probe of Ps chemistry, the positron delivers information about 
processes near and inside its terminal blob, while radiation chemistry primarily investigates the 
processes in isolated spurs. The main difference between the spur and blob comes from the 
factor of 10 difference in the initial number of ion-electron pairs they contain. 

In radiation chemistry, Onsager's theory of the geminate recombination is adequate for 
interpreting free-ion yields and their field dependence. On the contrary, the blob model is most 
appropriate for considering processes, in particular, Ps formation in the blob. At a limiting 
case, the theory we have developed gives a relation which only formally resembles the well- 
known Onsager equation of geminate recombination. We generalized the prescribed diffusion 
method for consideration of the processes in the presence of an external electric field. 

Thermalization of knocked-out electrons differs from the slowing of the positron. Being 
affected by an electric field of the parent ion, the electron becomes thermalized at a shorter 
distance than the positron, which easily escapes from its electrically neutral center of the blob. 
Thus ap becomes larger than the initial ion-electron separation, tq, in a spur. This is also 
related to the higher mobility of the thermalized positron compared to the electron mobility. 

An investigation of Ps formation in the presence of an external electric field enables us to 
better understand the behavior of intrablob electrons and the positron - peculiarities of their 
thermalization, quasi-neutrality condition, positron out-diffusion from the blob, formation of 
the qTiasifree positronium and its transformation to the bubble state. 

One of the interesting findings in our study is that the competition of the Ps formation with 
the other processes (like annihilation of positrons in the free state and lER) leads to a zero 
derivative of Pps{D) ai, D — 0. Another surprising but not very unexpected observation is that 
lER cannot compete significantly with the Ps formation. We found that the ratio of the Ps 
formation rate constant to the lER rate constant fcep/^ie is about 10^-10^. Different conditions 
of absorption of released energy in lER and Ps formation in the e"*" blob may be responsible for 
this result. 

In summary, we stress that experiments on Ps formation in the presence of an external 
electric field are extremely informative both for Ps chemistry and for radiation chemistry, and 
should be extended to higher fields and other substances. 
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APPENDIX A: ENERGY DEPOSITION AND TRACK STRUCTURE OF FAST 

POSITRON 

The average initial energy of fast positrons emitted from radioactive nuclei ^^Na or ^^Cu, fre- 
quently used as a source of positrons, is several hundreds of keV. Moving in a liquid, the positron 
loses about half of its kinetic energy in rare head-on collisions, knocking out 5-electrons. Tracks 
of these electrons form "branches" around the positron trail. The other half of energy is spent 
in numerous glancing collisions with molecules. Average energy loss in such a collision is 30-50 
cV (at maximum 100 eV). A secondary electron knocked out in a glancing collision produces, 
by turns, a few ion-electron pairs inside a spherical microvolume, called a spur in radiation 
chemistry. Its radius, asp, is determined by thermalization of the knocked out electrons in the 
presence of the Coulombic attraction of parent ions. Strictly speaking, a^p is the dispersion of 
the Gaussian distribution function, which corresponds to the end of thermalization process of 
the knocked-out e~. Based on^^''^^, the most probable value of Ugp in water is estimated to be 
f« 30 

While positron energy W > Wcyi, mean distance between adjacent events of ionization 
produced by the positron is greater than spur size 2asp (Fig.2).^° This means that spurs are 
separated from each other. The motion of a high energy positron is a quasi-straight line because 
li is less than the positron transport path Itr, which is the mean distance passed by the positron 
before it changes the initial direction of its motion by 90° (detailed in®''). When li < 2asp < kr 
or Wbi < W < Wcyi spurs overlap, forming something like a cylindrical ionization column. At 
the end of the track, when e+ energy becomes less than blob formation energy Wu (~ 500 
eV), 2asp becomes the largest parameter: 2asp > k, kr- This means that the positron starts to 
create a blob containing a few tens of ion-electron pairs because the average energy required 
for the formation of one ion-electron pair is 16-22 cV.^ Diffusion motion of the positron in 
the blob becomes more pronounced with decreasing energy. The positron frequently changes 
its momentum due to elastic scattering and the ionization of surrounding molecules. Roughly 
speaking all intrablob ionizations are confined within the sphere of the radius au (detailed 
definition sec below). The positron finally becomes subionizing and therefore its energy loss 
rate drops down by almost 2 orders of magnitude."'^ 

To gain some insight into typical values of the parameters involved we shall start with 
estimations of k{W) and ltr{W). Calculation of k{W) = Wiep/LET{W) is based on the data 
on LET of e='= (see, for example. Fig. 3. 19 in^^ or Fig. 5 in^*^). Usually the right side of the 
Bragg peak is well described by the Bethc formula for ionization slowing down, but its left (low 
energy) side strongly depends on corrections to the Bethe equation which consist in truncation 
of the dipole oscillator distribution at the maximum transferable energy.®^ From the relationship 
k{Wcyi) ~ 2asp we obtain that Wcyi ~ 3 keV, Fig. 3. 

Estimation of the transport path can be done in the frameworks of the Born approximation 
(wavelength of c"*" with the energy > 100 eV is small in comparison with the size of molecules). 
Below 1 keV /fr(M^) is mainly determined by small- angle electron-positron elastic scattering 
(the Rutherford part of the cross-section, related to e"*" scattering on nuclei, becomes important 
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above 1 keV):^^ 

ltr{W) = \—, atr{W)= f \fB{e)\\l-cose)27rsmede. (Al) 

The Born amplitude /b is calculated simulating a molecule of the liquid by an equi-electronic 
hydrogen-like atom^'^°. Energy dependence of Itr is shown in Fig. 3. 

At low energies positron scattering becomes more and more efficient and we must regard 
e"*" motion as diffusion-like. Itr is then considered as the energy-dependent mean free path 
between successive "collisions", which completely randomize the direction of the velocity of 
the particle. If the probability to pass distance r without such a collision is exp(— r/Zj^), the 
average squared distance is 2/^,.. After n collisions, the mean square displacement is nr^. 
Calculation of the same quantity assuming that particle propagation is governed by the usual 
diffusion equation gives nr^ = 6Dpt. Thus, we obtain the diffusion coefficient of the positron as 
Dp{W) = I'lr/'^'^ = hrVp/'i, where Vp is c+ velocity and r = t/n = f/vp = kr/vp is the average 
time between subsequent collisions. Integrating the relationship 

d{r^) = 6Dpdt = 2ltrVpdt = 2ltrdx = 2^^^— — (^2) 

\(lvV / (IX I ion 

from the energy Wi down to Wf, we obtain the diffusion displacement Rsd of e"*" during its 
ionization slowing down within this energy interval Wf <W <Wi. 




Now we are ready to define the blob formation energy Wu and the "radius" au of the blob. 
These quantities are determined from the following equations: 

ltr{Wu) = au, au = RsdiWu, Ry) - au- {AA) 

Their physical meaning is clear from Fig.4. Terminal blob is a spherical microvolume which 
confines the end part of the positron trajectory, where ionization slowing down is the most 
efficient (thermalization stage of subionizing positron is not included here). Mathematical for- 
mulation of this statement is twofold. Just after the first blob formation "step" (the thick arrow 
in Fig.4), ItriW^)) e"*" reaches the center of the blob. After that, slowing down displacement of 
the positron, Rsfi{Wbi,^j) — a^ should be equal to the "radius" of the blob, a^i, i.e. the blob 
in average embraces exactly the end part of the e"*" ionization slowing down trajectory. 

Solution of these equations is unique and shown in Fig. 3 for the case of liquid water. It is 
seen that numbers Wu ~ 500 eV and au = 40 A fulfill Eqs.(A4). Actually, we must proceed 
with all such calculations for each particular liquid investigated by the positron spectroscopy. 
However, one may assume that the values of au, Wu and other parameters related to slowing 
by ionization do not differ significantly from one liquid to another, because variations in the 
ionization potential and average electron density are small (Table I). 

At the end of slowing down by ionization and electronic excitation, the spatial distribution of 
the subionizing positron coincides with the distribution of the blob species (i.e. ~ exp(— r^/a^J). 
Further, during thermalization, e+ distribution becomes broader (about 30 A in water and 100- 
200 A in hydrocarbons), and its total dispersion is expressed as follows: 

/ dW \ 
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Estimation of requires knowledge of the stopping power, \dW/dx\mh, of a given liquid to- 
wards excitation of vibrations, scattering and energy loss properties of subionizing e^, and the 
spectrum of its initial energies Wq after the last ionization event. {. . .)wo denotes the average 
over Wo.^°'^^ Thus, contrary to the parameters related to ionization slowing down, Up strongly 
depends on the properties of each particular liquid. 



APPENDIX B: MICROSCOPIC QUASI-EQUILIBRIUM CONDITION FOR IONS 

AND ELECTRONS IN THE BLOB 

Let us estimate how much the distributions of electrons and ions differ at quasi-neutrality 
(quasi-equilibrium) condition which is achieved immediately after the electrons adjust them- 
selves to the current distribution of ions. As discussed before Eq.(3), it is reasonable to use 
Gaussian distribution functions to describe the spatial distribution of ions and electrons by the 
end of the thermalization stage: 

n„G.(r) = n ^^lV^ noG.(r) = no^22tl!M. (fii) 

no is the number of ion-electron pairs in the blob. We expect that, due to strong electrostatic 
attraction, the difference Aa = ae — ai is very small in comparison with the blob radius au ~ 
~ ag. The distribution of electrons is slightly broader than that of ions because electrons are 
much more mobile. 

Integrating the equation on electric field E(r) 

divE(r) = ^^(G,-Ge), {B2) 
e 

we find the electric field, produced by electrons and ions in a liquid: 

N 2 enor ^l"'^" Aa^au 4 euo Aa r f r''\ 
t][r) = — -- ^/xeiq)[-x)dx — > — 5- exp — 2" ■ (-^3) 

At the quasi-equilibrium, the electron flux should equal to zero everywhere. However, here we 
use the approximate (gaussian) functions instead of true distribution functions, so we can fulfill 
this condition only for a given value of r, for example, at r = a^'- 

CebeE{au) + DeWceiau) = 0, Ceiau) = noGeiau), be = (-B4) 

From this equation in the limit — <^ au, we find difference Oe — a^. Using Eq.(Bl) and 
Eq.(B3), we obtain 

Vceiau) = Ce{au)-, E{au) = % , e = 2.718 .... (55) 

Substituting these equations into Eq.(B4), we arrive at 

Aa c r. = ^. {B6) 

rioTc sT 

Numerically, it yields Aa ~ 1 A (an estimation of no and au is given in Appendix A). The 
potential, which corresponds to Eq.(B3), is 

ip{r) = A^^exp(-rVa^,) = — exp(l - rVa,^,), -Wip^r) = E(r). {B7) 
y/neabi au ere 
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Note that the apphcation of an external field of about 30 kV/cm to a liquid does not 
significantly perturb the ion-electron distribution in the blob. It leads to less than 1 A shift of 
the electron distribution with respect to the ions. 



APPENDIX C: AMBIPOLAR DIFFUSION EXPANSION OF THE BLOB 

Even neglecting ion-electron recombination, the out-diffusion of the blob species (electrons 
and ions) can not be considered independently because of the strong electrostatic interaction 
between charged particles. The intrablob electric field efficiently suppresses out-diffusion of the 
electrons but increases (duplicates) the ion diffusion coefficient.®^ Coupled diffusion equations 
are written as 



dci{r,t) 
dt 

9ce(r, t) 



= A (Aq + ^Av^) , (Cla) 

= De (ACe - ^A<^) . (Clb) 

Di and are the local diffusion coefficients of ions and electrons and A stands for the Laplace 
operator. Our aim now is to relate them to the kinetics of blob expansion. As before, Eqs.(3) 
are used as initial conditions. To Eqs.(Cl), we add the Poisson equation for an electrostatic 
potential (p{r): 

A^=-^(Q-Ce). (C2) 

In Appendix A, we demonstrate that the spatial distribution of electrons is very close to that 
of ions, so 5c = Cj — Ce -C Cj or c^. Using Eq.(C2), we evaluate the second terms in the right 
hand side of Eq.(Cl) as follows: 

—Aip ~ incrcdc, c(r) = Ce ^ Ci, —. 

Ace and Acj are on the order of c/a^^. Their difference, Ace — Acj ~ Sc/ali, is negligible 
compared to the second term, ATrcrcSc, on the right hand side of Eqs.(Cl), because their ratio 
is small: 

6c/ all _ abi ^ ^ 



We assumed c ~ no/a^^. 

Moving fast {D^ » Dj), electrons rapidly adjust themselves to the current distribution of 
ions. So dc^jdt — 0, which gives 

^A(^ = Ace. (C3) 



Together with Ace ~ Acj and ^A(/7 ^A(/?, we arrive at 

-A(^ = Aci. 



T 

QC 

~T 

This means that the term Acj is duplicated in Eq.(Cla). The motion of ions in the blob thus 
obeys the simple diffusion equation, but with a twice larger diffusion coefficient: 
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^ = (C4) 

This process is called ambipolar diffusion and Damb — '^Di is the ambipolar diffusion coefficient. 
Half of it is related to the diffusion of ions and the rest is due to the electric field of blob electrons. 
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Figure Captions 



Figure 1 

Dependence of o-Ps intensity Is — 3Pps/4- 100% on external electric field D in different liquids: 

□ - cyclohcxanc {Ce,B.i2), © - n-hexane (C6H14), A - isooctanc (Cgliis)- ^ ^ hexafluorobenzene 
(CeFg) and O - benzene (CqH^). The size of symbols reflects a statistical uncertainty only. 
Broken lines represent the fit according to Eq.(22), which does not account for the free positron 
annihilation. They solely show the hmiting case of the theory. Sohd hues show the fit when 
positron annihilation with the hfetime T2 and lER {Wie/Wep — 0.003) are taken into account 
in Eq.(32). The lines which account only effect of the positron annihilation lie in between 
respective dashed and solid lines and not shown here. Experimental data suggests the presence 
of lER on the level Wie/Wgp < 0.005 or positron localization with a rate lower than I/T2. 
Corresponding parameter values are shown in Table II. 

Figure 2 

Scheme of the end part of the e+* track and Ps formation. 
Figure 3 

Dependence of the mean distance between adjacent ionizations li{W) and the transport mean 
free path ItriW) of the positron on its kinetic energy."^° Spurs start to overlap when li{Wcyi) ^ 
2asp, forming cylindrical-like column (the arrow at right). Formation of the blob begins when 
the positron energy fulfills equation (A4): ltr{Wbi) = Rsi{Wbi, Ry)/'2. 

Figure 4 

Schematic view of the terminal positron blob. Positron motion is simulated as random walks 
with the energy dependent step, ltr{W)- 
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TABLE I. Parameters of the liquids 
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i?ws is the radius of the Wigner-Seitz sphere at room temperature, AttR^^/S = 1/n, where n 

is the molecular concentration. 

£ = £00 is the dielectric permittivity of the liquid.™ 

Z is the total number of electrons in molecule 

He = Zn: Average electron density of liquid 

T3: Ortho-Ps lifetime.^^ 

T2'- Free-e^ lifetime. 

be'. Mobihty of an excess electron.^^'^^ 

bp-. Positron mobility. In^^ e"*" mobility was defined as a coefficient of proportionality between 

e"*" drift velocity and external electric field D. Thus, in our notation respective numbers should 

be multiplied by a factor of e. 

To'- Ion-electron initial separation in a spur.*^^'^^ 

Vf^: Electron work function. 

Gas- phase ionization potential (/g)-^° 

Liquid phase ionization potential (II)-'^^ 

Reference^^ 
'^^ Reference^ 
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TABLE II. Parameters obtained from the fit of the data 
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3.1 • 10^^ 
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CgFg 


68(4 

V / 


0.687(3) 
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5.3 • 10^3 






79(4) 


0.709 3 
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6.3 • 10^3 
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0.0011 
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185(4) 
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0.0006 


8.3 • 10^^ 


7.1 • 10^° 


cyclohexane 


196(4) 


0.326(3) 




7.7 • 10^^ 
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0.336(3) 


0.0008 


8.5 • lO^'' 






222(5) 
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0.0018 


9.1 • lO^'' 


- 




226(5) 


0.344(3) 


0.0010 


9.3 • 10^^ 


4.4 • 10^° 


benzene 


166(3) 


0.410(3) 




1.5 • 10^^ 




CeHg 


214(5) 


0.451(3) 


0.0043 


2.1 • 10^^ 
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0.484(3) 


0.0112 


2.5 • lO^'' 






224(5) 


0.461(3) 
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2.2 • lO^'' 
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The first line for each hquid corresponds to the fit based on Eq.(22) (dashed hnes in Fig.l), which 
does not account for free positron annihilation. We list these numbers to show the simplest 
limiting case of the theory. These numbers differ from corresponding data in^^ because the 
field-dependent bias of subionizing positron during its thermalization is taken into account in 
the present work. 

Numbers on the second line show the fit including positron annihilation with lifetime T2 based 
on Eqs. (32-33), but lER was not taken into account. 

Data on the third line corresponds to the same case as above, but positron localization is 
assumed with the rate l/Tp°'^ = I/T2. Respective curves for Pps(-D) in the last two cases are not 
shown in Fig.l. They are in between the dashed and the solid hnes. 

Numbers on the last line of each group are obtained based on Eq.(32), where, together with the 
positron annihilation with the rate T2^, a fixed fraction of lER is introduced: Wie/Wgp — 0.003 
(maximum possible level still acceptable by experimental data). 

/ugp and /c^g are recalculated from Jal^ + yj/^^ and Wie via Eq.(21) and Eq.(31). It is assumed 

that no = 30, = 40 A and Dp = Tbp/e, where bp is taken from Table I (e"*" mobility for 
cyclohexane is assumed to be the same as for hexane). 
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